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ABSTRACT 

Based on detailed analysis of radio and X-ray observations of a flare on 2002 April 11 aug- 
mented by realistic 3D modeling, we have identified a radio emission component produced 
directly at the flare acceleration region. This acceleration region radio component has dis- 
tinctly different (i) spectrum, (ii) light curves, (iii) spatial location, and, thus, (iv) physical 
parameters from those of the separately identified, trapped or precipitating electron compo- 
nents. To derive evolution of physical parameters of the radio sources we apply forward fitting 
of the radio spectrum time sequence with the gyrosynchrotron source function with 5 to 6 
free parameters. At the stage when the contribution from the acceleration region dominates 
the radio spectrum, the X-ray- and radio-derived electron energy spectral indices agree well 
with each other. During this time the maximum energy of the accelerated electron spectrum 
displays a monotonic increase with time from ~ 300 keV to ~ 2 MeV over roughly one minute 
duration indicative of an acceleration process in the form of growth of the power-law tail; the 
fast electron residence time in the acceleration region is about 2 — 4 s, which is much longer 
than the time of flight and so requires a strong diffusion mode there to inhibit free-streaming 
propagation. The acceleration region has a relatively strong magnetic field, B ^ 120 G, and 
a low thermal density, <2 ■ 10^ cm~^. These acceleration region properties are consistent 
with a stochastic acceleration mechanism. 

Subject headings: Sun: flares — acceleration of particles — turbulence — diffusion — Sun: magnetic 
fields — Sun: radio radiation 



; 1. INTRODUCTION 

Probing the acceleration region of solar flares is 
\ known to be a highly nontrivial tas k ([Bastian et al.l 
. [200llXu et"all[200l iHolmanllMl since the direct 
emission from the acceleration site is often weaker 
than other competing emissions, e. g., soft X-rays 
(SXRs) from thermal plasma and microwave con- 
tinuum from a magnetically trapped fast electron 
\ component. Rece ntly, we identified and reported a 
cold, tenuous fiare (jFleishman et al.ir2011| ) , which dis- 
played neither hot coronal plasma nor a magnetically 
I trapped population of fast electrons in a coronal loop. 
This rare but favorable combination of fiare proper- 
ties allowed us, for the first time, to firmly identify, 
spectrally and spatially, the radio emission compo- 
nent produced directly at the acceleration region and 
then derive physical parameters of the acceleration 
site and accelerated electron population. 

Although detection of even a single acceleration 
region is important for improving our understand- 
ing of where and how the flare electrons are acceler- 
ated, it might seem discouraging because such 'clean' 
cases without plasma heating and electron trapping 
are extremely rare. In contrast, in a typical case, 
the flare plasma heating is an essential component, 
vividly seen in coronal SXR sources, while in the mi- 

^ Center For Solar- Terrestrial Research, New Jersey Insti- 
tute of Technology, Newark, NJ 07102 

^ loffe Physico- Technical Institute, St. Petersburg 194021, 
Russia 

^ Department of Physics and Astronomy, University of 
Glasgow, G12 8QQ, United Kingdom 



crowave range three source contributions (separate 
regi ons of acceleration, trapping, and precipitation, 
e.g.. IAschwandenlll998[ ) may be present. To study 
acceleration processes, it is necessary to distinguish 
the acceleration region contribution in the presence 
of these other two competing emissions. An ideal 
way to make this distinction, at least for cases when 
the acceleration is spatially displaced from the other 
components, is through the use of high spatial and 
spectral resolution microwave imaging spectroscopy, 
which is not yet routinely available. In the mean- 
time, however, we are restricted to selecting rare cases 
where the direct emission from the acceleration region 
can be separated spectrally and/or temporally from 
the competing contribution in the spatially integrated 
dynamic spectrum. 
Theoretical consideration of radio emission ex- 

J iected from the acceleration regions in flares 
Li fc Fleishi^[2009t iPark fc FleishmanllMol ) sug- 
gests that the corresponding contribution should peak 
at a few GHz and have a reasonably narrow spec- 
trum. This is, in particular, the case of the radio 
spectrum fro m the acceleration reg ion of the cold, 
tenuous flare (jFleishman et al.ll2011 ): however, hav- 
ing a spectral peak at this decimeter range does not 
guarantee that the emission originates from the accel- 
eration region, so some additional evidence is needed. 
In this paper we present the study of a flare in which 
a certain fraction of the radio emission (at a some- 
what narrow spectral range during a limited time) 
can be confldently attributed to the acceleration re- 
gion, even though the competing contributions are 
strong and overall comparable with or even dominat- 
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Fig. 1. — Overview of April 11, 2002 flare: Phocnix-2 and 
OVSA dynamic spectra, top panels. GOES (3 s) lightcurves 
as measured by GOES-10 spacecraft. RHESSI (2 second bins) 
lightcurves (bottom panel) in: 3-9 koV (black), 9-15 keV (red), 
15-30 keV (orange), 30-100 keV (blue). 



ing that from the acceleration region. In particular, 
we will demonstrate that the radio detection of the 
acceleration region in the given event is favored by 
stronger magnetic field at the acceleration site com- 
pared with the coronal 'trapping site', where the fast 
electron accumulation occurs at a later stage due to 
the well-known effect of magnetic trapping. Thus, the 
gyrosynchrotron (GS) emission from the acceleration 
site dominates until the number of the magnetically 
trapped electrons rises above the level needed to dom- 
inate over the acceleration region contribution. 

In section 2 we present the observations from 
Reuven Ramaty H igh Energy Sola r Spectroscopic 
Imager (RHESSI, iLin et alJ [200l in the X-ray 
range, Owens Valley Solar Array (OVSA) at 1- 
18 GHz, and Phoenix - 2 spectrometer at 0.1-4 GHz 
(|Messmer et al.l I1999D . Context SXR observations 
made by GOES-10 and SOHO/MDI measurements 
of the photosphcric magnetic field are utilized as well 
along with a new 3D modeling tool, GX_Simulator, 
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Fig. 2. — RHESSI spatially integrated X-ray spectrum above 
6 keV for two HXR peaks of the event: 16 : 19 : 43 - 16 : 
22 : 01 UT - top panel and 16 : 23 : 07 - 16 : 25 : 13 UT - 
lower panel. In both panels, brown line indicates background, 
red-dashed lines thermal and non-thermal components, blue- 
line is the total fit (thermal plus non-thermal), and black and 
brown histograms X-ray data for the signal and the background 
respectively. 

which we have developed and recently included into 
Solar Softwarfl In section 3 we present spectral fit- 
ting of the radio and hard X-ray spectra using 5 to 
6 physical parameters, and show that the accelera- 
tion region can be separately identified. We outline a 
possible 3D geometry in section 4 and conclude with 
a discussion of the results in terms of a stochastic 
acceleration mechanism. 

2. OBSERVATIONS 

2.1. X-ray Imaging and Spectroscopy 

The flare demonstrates two hard X-ray peaks as evi- 
dent from Figure [H The spatially integrated RHESSI 
X-ray spectrum (Figure [2]) over the duration of each 
peak indicates the presence of both thermal and non- 
thermal components as often seen in X-ray spectra 
(e.g. see Kontar et al. 2011, as a review). To achieve 
better count rate statistics and energy resolution, we 
summed all front segment detectors but the detec- 
tors 2 and 7 were avoided (see Smith et al. 2002 for 
the d etails). Spectral ana lysis was done using OS- 
PEX (jSchwartz et al.ll2002D with systematic errors set 
to 0.02%. As the flare appeared on the solar disk 
{x ~ 500", y ~ - 180" ; heliocentric angle ~ 56° 
), albedo correction (jKontar et al.ll2006h was applied 
assuming isotropic emission (minimum correction). 
The spectrum was fitted with a standard thermal plus 
non-thermal thick-target model (Figure [5]), wi th the 
i sotro pic bremsstrahlung cross-section following lHaud 
(11993). 

The temporal variation of plasma parameters ob- 
tained from the HXR fit with time is presented in 
Figure [31 The derived accelerated electron spectral 

|ww. Imsal ■ com/solarsof t/ssw_packages_lnf o .html | 
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Fig. 3. — Spectral fit parameters from the thermal plus thick- 
target X-ray model for 20 second time intervals of RHESSI data 
(black curves) as well as emission measure and temperature as 
obtained from GOES-10 data (blue curves). 



evolution demonstrates typical soft-liard-soft behav- 
ior and shows two clear peaks in hard X-rays, well vis- 
ible in 30 — 100 keV lightcurves (Figure [T|). The emis- 
sion measure of the thermal plasma increases from 
< 10''^'' cm~^ to ~ 10*® cm~'^, while the temper- 
ature of the plasma stays rather constant, ranging 
from ~ 20 to ~ 15 MK over the course of the flare. 
The emission measure and the temperature observed 
by GOES-10 are somewhat higher and lower, in the 
ranges above 2-10*® cm~'^ and between 12 and 10 MK, 
respectively, but show a similar temporal behavior. 

RHESSI imaging (Hurford et al. 2002) reveals a 
spatially well-defined soft X-ray (6-12 keV and 12-20 
keV) source (Figure 2]) , visible over the entire time 
of the flare, whose location gradually moves east- 
ward by roughly 10" over the flare duration. HXR 
sources imaged in the 20-40 keV range have higher 
temporal variability, which coincides spatially with 
the SXR source at the troughs of X-ray light curves 
at 16 : 19 UT and 16 : 23 UT but deviates from 
this location at the peaks, when the electron ac- 
celeration rate is the highest. Based on SXR-to- 
HXR spatial relationship alone, one could interpret 
the event in terms of the standard flare picture, i.e. 
a looptop SXR source and two uneven HXR foot- 
points. This interpretation, however, is supported 
neither by the magnetic field data shown in Figure 
[5] nor by the observation-based 3D modeling given 
below in Fig. [lOl We argue based on analysis of 





Fig. 4. — Spatial distribution of X-ray emission from April 
11, 2002 flare in various energy ranges. IBackground image 6 — 
12 kcV, 12-20 keV green (contours are at 50, 70, 90 % levels), 
and 20-40 keV dark blue (contours are at 70, 80, 90 % levels). 
The X-ray images are reconstructed using clean algorithm. 
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Fig. .5.— 20-40 keV HXR contours [45; 60; 75; 90%] for 
16:20-16:21 UT superimposed on the SOHO/MDI line of sight 
magnetic map. It is clearly seen that the HXR source is pro- 
jected to the photosphere region free from any noticeable mag- 
netic field enhancement, which supports the idea that the HXR 
source is coronal, rather than a footpoint. 



the microwave data, which we will show later, that 
the main HXR source remains coronal most of the 
time, but support for this can be seen in the fact 
that the 20-40 keV source in Fig. [5] does not over- 
lie any concentration of magnetic field, that could 
correspond t o foot points often seen in flares (e.g. 
Kosud et afl 119921: lEmsli e et al.l [20031 iKontar et atl 
20101 : iBattagha fc K ontar 2011[). The only time pe- 
riod when the main HXR source possibly could be 
identified with a footpoint is late in the first peak 
16:21- 16:22 UT when the 20-40 keV source is more 
strongly displaced to the north from the SXR source, 
while later the HXR source moves back slightly to the 
south. A weaker HXR source, seen as a south-east ex- 
tension of the main HXR source in Figures [S] and [HI 
left, could correspond to a footpoint; alternatively, it 
can correspond to another part of the same magnetic 
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Fig. 6. — Left: Evolution of the spatial brightness peak 
of the radio emission at 2.6-4.2 GHz from April 11, 2002 
flare. Background: HXR image at 20 - 40 keV at 16:20:00- 
16:21:00 UT. Symbols are OVSA image centroid positions sep- 
arated by 8 s time interval (OVSA temporal resolution for 
imaging): green triangles are for snapshots from 16:19:55 to 
16:20:15 UT, pink asterisk for 16:20:20 UT, black triangles 
from 16:20:25 to 16:22:00 UT, white triangles from 16:23:40 
to 16:24:20 UT, and the blue asterisk is for a late decay phase 
of 16:25:08 UT; larger triangles correspond to later snapshots 
within each group. The sequence of the contours clearly in- 
dicates that the radio source is located at the northern part 
of the HXR image and stays there during the entire impulsive 
phase, then moves southward to exactly match the HXR cen- 
troid position and stays there the entire decay phase of the 
first peak, then it moves westward during the second peak, 
and finally returns back to HXR centroid location at the late 
decay phase. The accuracy with which the contours with the 
same color coincide with each other indicates the accuracy of 
the position measurements. The synthesized bean is shown by 
the dashed white oval. Right: centroids and 95% contours of 
the three dominant source (corresponding to the green, black, 
and white sources on the left panel) positions superimposed on 
the SXR image; blue dashed contour shows the HXR, source at 
60% level. 

structure or to a lower, more compact loop. 

2.2. OVSA imaging 

We employ OVSA frequency synthesis imaging 
around the spectral peak frequency, 2.6-4.2 GHz, 
while no reliable imaging was possible at the higher 
frequencies, ^ 10 GHz due to the low flux levels. 
Given a limited OVSA spatial resolution at low fre- 
quencies, the radio sources were unresolved, although 
useful information on the source location at different 
stages of the burst has been obtained. 

Figure [5] displays the centroid location of the radio 
emission during different stages of the event superim- 
posed on the HXR and SXR images plotted for time 
interval 16:20:00-16:21:00 UT. A few important con- 
clusions can be made based on the presented spatial 
relationships. First, we note a clear spatial evolution 
of the radio source, which is initially (during the im- 
pulsive peak of low-frequency microwave emission at 
2.6-4.2 GHz) displaced by roughly 10" from the ac- 
companying HXR source with even greater displace- 
ment, by roughly 20", from the accompanying SXR 
source. Upon transition to the decay phase, the ra- 
dio brightness peak shifts by ~ 10" to exactly match 
the brightness center of the HXR source, which re- 
mains displaced 10" from the SXR source. Finally, 
during the second temporal peak of the flare around 
16:24 UT, the radio brightness peak is displaced by 
10" to the east, after which it eventually returns back 
to the location of the radio source during the decay 
phase of the first flare peak. 

2.3. Radio to X-ray Timing 
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Fig. 7. — Radio to HXR timing with Phoenix-2 and OVSA at 
different frequencies. Black thin curve, same throughout all the 
panels: HXR (20-40 keV) light curve with 2 s resolution nor- 
malized to 1; blue thick curves: resampled (to the RHESSI 2 s 
resolution) Phoenix-2 light curves at different frequencies in- 
dicated at the panels; green curves: spline interpolated OVSA 
light curves from original 4 s resolution to the RHESSI 2 s res- 
olution at different frequencies indicated at the panels. The 
insets: corresponding Radio-to-HXR lag correlations in which 
the time delay is printed by either blue (for Phoenix-2 data) or 
green (for OVSA dat a ). T he light curves at other frequencies 
are given in Figs. [TTI— 1151 in the end of the file. 

Study of the emission light curves including cross- 
correlations and delays between various microwave 
and X-ray channels is often highly helpful in iden- 
tifying distinct emission components, which as noted 
in the introduction may be produced by electrons at 
the acceleration site, by magnetically trapped elec- 
trons, or by precipitating electrons. Figure [7] displays 
temporal relationships between radio light curves 
recorded by Phoenix-2 and OVSA on one hand and 
HXR light curves recorded by RHESSI on the other 
hand. To perform the lag cross-correlation analy- 
sis between the radio and HXR light curves, the 
Phoenix-2 light curves were resampled, while OVSA 
light curves spline interpolated to match 2 s RHESSI 
resolution. The displayed radio light curves are dis- 
tinctly different from each other depending on the 
frequency and can be grouped as follows. 

The low frequency radio light curves, / < 2 GHz, 
display a significant time delay compared with the 
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HXR light curve: up to half a minute as determined 
from lag cross-correlation and up to 45 s if the peak 
times are considered. Then, at higher frequencies, 
2 GHz < / < 5 GHz, there is no significant delay be- 
tween the light curves; the delay varies within ±2 s in 
this frequency range. However, the various radio light 
curves do show spectral evolution within this range: 
they become more and more impulsive as frequency 
changes from 2 to 2.8 GHz, then the light curves stay 
roughly similar to each other at 2.8 — 3.8 GHz, and 
then they again become less impulsive at 4 — 5 GHz. 
At higher frequencies, 5.0 — 7.5 GHz, the radio light 
curves are again significantly delayed compared with 
the HXR light curve similarly to the lowest frequency 
range. Finally, the highest frequency range is char- 
acterized by light curves with an impulsive peak co- 
inciding with the HXR peak and with the impulsive 
radio peak at frequencies 2.8 — 3.8 GHz, but followed 
by a more extended part of the light curves. 

We conclude that the radio light curves consist 
of two essentially different components — impulsive 
and delayed — whose relative contributions depend 
strongly on the frequency0 The delayed radio 
light curves behave exactly as expected for the mi- 
crowave emission produced by a mag netically trapped 
component of fast e lectrons (e.g. , iMelnikovl 119941: 
Bastian et all Il998t iKundu et all 120011: IMelnikovl 
2006( ). In contrast, the impulsive radio light curves 
around the spectral peak frequency of 3.2 GHz reach 
maximum flux simultaneously with the HXR light 
curve and so cannot be ascribed to the trapped 
component. Instead, their timing shows remark- 
able similarity to that foun d in the cold, tenuous 
flare (jFleishman et al.ll20H"l ). I.e. microwave emis- 
sion produced directly from an acceleration region 
at low frequencies and from precipitating electrons 
at high frequencies. We suggest that the impulsive 
and delayed emissions come from physically different 
sources, which is consistent with their different spa- 
tial locations discussed in the previous section. 

3. RADIO SPECTRAL FIT 

An alternative way of looking at spectral evolu- 
tion, rather than frequency- dependent light curves, 
is through time- dependent spectra during the burst. 
These are shown in Figure HI which includes both the 
measured total power radio spectra (points) and spec- 
tral fits to be described below. The observed spec- 
tral evolution is in marked contrast with the even t 
of 2002 July 30 analyzed by Flcishni^eEaD (|2011| ). 
where almost no spect r al evo lution was present. In- 
deed, iFleishman et al.l ()2011[ ) successfully fitted the 
entire sequence of the radio spectra in the 2002 July 
30 event with a source model having only one free 
parameter — the instantaneous number of the acceler- 
ated electrons at the source. Judging from the sig- 
nificant spectral evolution in the event under study, 
such a simplified approach to the radio spectral fit 
cannot work. A multi-paramet e r fit s imilar to that 
developed by IFleishman et al.l ()2009() is needed in 

^ No impulsive component (acceleration region contribution) 
is detected during the second peak of the flare occurring around 
16:24:00 UT; not shown in Figure [7] 
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Fig. 8. — Radio spectral fit as described in the text (solid 
curves) of the OVSA spectra (symbols with scatter) at conse- 
quent snapshots with 4s cad ence; Part I. Other spectra with 
fits are given in Figs. [TBI— 1191 below. 

this case. However, IFleishman et al.l (|2009D applied 
their multi-parameter fit to a sequence of spatially 
resolved (model) spectra, while here we are forced to 
deal with the total power spectrum, in which source 
nonuniformity can be e xpected to play a major role 
(|Kuznetsov et al.ll201lD . 

Note that many of the spectra in Figure [5] show 
two spectral peaks. A similar two-component struc- 
ture was present in all radio spe ctra in the 2002 July 
30 event (|Fleishman et al.ll201lD . which was success- 
fully interpreted as a combination of a main (low- 
frequency) coronal source, representing in fact the 
acceleration region in that event, and a secondary 
(higher-frequency) source formed by the precipitat- 
ing component of the fast electrons. Here we employ 
a similar two-component spectral model, i.e., that the 
radio spectrum is formed in two uniform sources phys- 
ically related to each other. While many of the radio 
spectra do not show two clearly separated spectral 
peaks, nevertheless, we fit two homogeneous spectral 
components at all times because they are needed to 
account for the source inhomogeneity, which broad- 
ens the spectrum to the extent that it cannot be fitted 
by a single-component, homogeneous GS spectrum. 

Specifically, we adopt the premise that the bulk of 
low-frequency emission comes from a coronal source 
with unknown (free) parameters such as magnetic 
field B\ and thermal density rie, where a power- law 
population of fast electrons is present with an un- 
known number density n^, spectral index 5, and high- 
energy cut-off i?max, but with a fixed low-energy cut- 
off -Emin = 17 keV, as shown in the lower panel of Fig- 
ure |3l The GS emission, especially its polarization, 
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depends on the viewing angl e and angular distribu- 
tion o f the fast electrons (e.g.. lFleishnian fc Melnikovl 
Since no polarization data are available, we fix 
the viewing angle at a quasi transverse value, 6 = 80° 
and adopt the simplest case of an isotropic electron 
distribution. 

Based on the imaging data we fix the main source 
size to be 25" in each dimension so the low-frequency 
source is fully characterized by these five free parame- 
ters. For the second ('precipitating') source we begin 
with one free parameter, the magnetic field B2, keep- 
ing 7ie, -Emin, and E'max valucs the same as for the 
'coronal' source. The spectral index of the instanta- 
neous distribution of precipitating electrons is taken 
to differ by 1/2 from that of the main electron popu- 
lation, Spr = (5 + 1/2, to account for the time of flight 
effect, TTOF L/v oc which assumes energy- 

independent escape and a free propagation transport 
regime for the precipitating electrons, and the size 
of the high-frequency source is corrected to conserve 
magnetic flux along the flaring magnetic flux tube. 
Certainly the adopted precipitation mode l may be 
oversimplified (see, e.g. Ivan den QordlfToonl : ISu et al.l 
1 2 1 It lHolmanl[2012i ) but it is sufficient for the spectral 
fitting and allows clear isolation of the low-frequency 
spectral component during the time when the double 
peak structure is clearly present in the radio spec- 
trum. Since our focus is on the low-frequency com- 
ponent, we forego more-detailed discussion of the pre- 
cipitating source. The gyrosynchrotron source func- 
tion is computed by the fast G S code developed by 
[Fleishman fc Kuznetsovl (I2010D . 

Although this six-parameter fitting succeeds 
throughout the entire flare duration, there is a mis- 
match between the fit and observed spectrum at some 
time frames at the lowest frequencies. Th is is a well- 
known indication of source nonuniformitv (iNita et al.l 
[200I lAltvntsev etldl [20i08t iQiu et al.ll2009[ ). which 
could have been accounted by adding a contribution 
from the electron accumulation site (see below); we 
avoid this complication, however, because the param- 
eters of this third spectral component could not be re- 
liably constrained from the fit. The derived physical 
parameters show noticeable scatter, particularly in 
the magnetic field B2, whose values clustered around 
several discrete levels near 800, 400, and 200 G dur- 
ing different fiare stagefl Accordingly, to improve 
stability of other derived parameters we fixed the 
B2 (whose exact value is inessential to what is dis- 
cussed below, but needed to obtain a good fit, es- 
pecially, for the spectra with two peaks) to one of 
the above values based on the original six-parameter 
fit, for time ranges corresponding to peak 1 impul- 
sive phase (800 G), peak 1 decay phase (400 G), peak 
2 rise phase (200 G) and peak 2 decay phase (again 
800 G), as is explicitly shown in Figure IHl We then 
repeated the fitting for the remaining five free param- 
eters in the low-frequency ('coronal') source, and the 
same parameters with the above-stated adjustments 
in d and source size for the high-frequency ('precipi- 

^ In fact, when B2 ~ 200 G, the contribution from this 
source is ncgUgible, thus, a one-component GS fit works well 
during this period of time. 
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Fig. 9. — Radio source parameters as derived from the OVSA 
spectral fit for five parameters of the low-frequency coronal 
source and adopted magnetic field value B2 for the 'precipi- 
tating' source as described in the text. A solid curve at the 
thermal plasma number density shows a number density evo- 
lution of the SXR source derived from emission measure from 
the RHESSI fit. Two top panels show the radio light curves 
recorded at 3.2 GHz and 9.4 GHz given for the reference pur- 
pose. 

tation') source, therefore, the presence of the precip- 
itation source does not add any additional free pa- 
rameter to the model. The time evolution of the de- 
rived parameters of the coronal source and adopted 
magnetic field B2 of the precipitation source is shown 
in the bottom six panels of Figure |9l compared with 
the observed 3.2 and 9.4 GHz microwave time profiles 
(upper panels). 

The derived evolution of the physical parameters 
deserves some discussion. The magnetic field in the 
low-frequency source is about 120 G during the impul- 
sive phase of the radio burst, while it drops quickly to 
30-50 G at the transition to the decay phase around 
16:20:20 UT. Remarkably, this magnetic field change 
derived from the spectral fit, Figure [9l happens at the 
very same time as the 10" shift of the spatial bright- 
ness peak, sec Figure |6l This implies that it makes 
sense to distinguish between these two spatially dis- 
tinct low-frequency sources, which as we show be- 
low represent the very acceleration region (the early 
source, with B 120 G, producing the impulsive ra- 
dio emission) and the classical looptop radio source 
(the later source, with i? ~ 40 G, producing the ra- 
dio emission from magnetically trapped electrons over 
the decay phase) spatially coinciding with the HXR 
source. 

Acceleration region. At the impulsive low- 
frequency source, - 16:20:00-16:20:20 UT, the ther- 
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mal number density obtained from the radio fit is 
somewhat low, Uf. < 2 ■ fO^ em~^, implying the ra- 
dio source is located in the corona, not at a chromo- 
spheric footpoint, while the number of nonthermal 
electrons is consistent with the acceleration rate de- 
rived from the HXR data, (1 — 3) • 10^'' electron/s, 
see Figure [31 if they reside at the radio source for 
2-4 s, which requires the strong diffusion transport 
mode. The radio derived electron spectral index does 
not display any significant departure from the HXR 
derived electron spectral index during this time inter- 
val, cf. Figure [3] All these properties are similar to 
those determi ned for the acceleration site in the 2002 
July 30 event (iFleishman et al.ll2011[ ). from which we 
conclude that we have here another instance of the 
acceleration region detection in a solar flare. The 
electrons accelerated at this source escape from there 
in roughly 3 s and then accumulate in another, 'trap- 
ping' source, which dominates the radio spectrum and 
spatial location after 16:20:20 UT. The acceleration, 
however, continues for a longer time: we note that the 
maximum electron energy, -Emax, displays a mono- 
tonic increase from ^ 300 keV to ^ 2 McV over this 
phase of the burst (~16:20:00-16:20:50 UT), which 
is reasonable to interpret as the growing of a power- 
law 'tail', i.e., the very process of the electron ac- 
celeration. Thus, the impulsive low-frequency source 
dominating the radio emission over roughly 16:20:00- 
16:20:20 UT, which produces fast electrons and sup- 
plies them to the coronal trapping site until at least 
16:20:50 UT, can confidently be identified with the 
acceleration region of the flare under study. 

Electron accumulation site. Transition to the 
gradual decay phascQ at about 16:20:20 UT manifests 
the stage when the trapping site has accumulated a 
sufficient number of fast electrons to dominate the 
radio spectrum. At this time the derived number 
of accelerated electrons with £' > 20 keV reaches a 
maximum of ~ 10^^ electrons, which corresponds to 
the number density of the accelerated electrons of 
rir ^ 2 ■ 10^ cm"'^ for the adopted source volume. 
For the determined earlier acceleration rate of 
(1 — 3) • 10'^^ electron/s, having a total electron 
number of ~ 10^^ requires a highly efficient electron 
trapping with the trapping time longer than 30 s. 
Indeed, that long trapping time is fully confirmed 
by the measured delay between the HXR/impulsive 
radio light curves and 'non-impulsive' radio light 
curves, see Figure 1111 some of which are delayed 
by almost one minute. This means that in this 
electron accumulation site the low-frequency radio 
emission is dominated by the magnetically trapped 
electron comp onent, whi c h is of te n seen in microwave 
burst s fe.g. iMelnikovl 119941: iMelnikov fc MagunI 
1991 [Banian et al.l 119981: iLeTfc Garvl 120001 



netic trapping implies a weak diffusion mode at the 
accumulation site, mediated, perhaps, by Coulomb 
collisions. 

Indee d, it is we ll known f e.g.. IMelrose k, Brownl 



19761: IM elnikov Iggi"; ''Bastian e t all Il998t 
Melnikov fc Magun.l998.:,Lee fc Garv.2000HMelmkml 



Kundu et all 120011: IMe lnikov et al.' '2002': 'BastianI 
2006t iTzatzakis et aLl lloOS: Reznikova et al. 2009;). 
In contrast to the strong diffusion mode in the 
acceleration region established above, effective mag- 

Its main parameters, Bi ~ 30—50 G and rie < 5-10^ cm"'', 
clearly indicate its coronal location, although spatially distinct 
from the acceleration region showing a different coronal loca- 
tion. 



20061 ) that Coulomb losses of the trapped electrons 
result in a flattening of the fast electron energy spec- 
trum. This effect must be present and accounted for 
since the Coulomb loss time is about 10 s for 20 keV 
electrons, which is much shorter that the trapping 
time at the decay phase of the burst. Indeed, the 
recovered evolution of the electron spectral index 
(Figure [His) displays a monotonic decrease over the 
entire decay phase, nicely in agreement with the 
theoretical expectation, while the maximum electron 
energy remains roughly constant at the level of 
~ 2 MeV implying that no acceleration to even 
higher energy is demanded by the radio spectrum fit. 
Although some of the details of the accumulation site 
interpretation are not unambiguously confirmed, e.g., 
some regimes of the wave scattering ca n result in a 
weak diffusion mode, (se e, e.g.. Bespalo v et al.lll987l : 
IStepanov fc TsanI I2002D . and still be compatible 
with efficient magnetic trapping, our main point 
is that this radio source (coinciding with the HXR 
source) is distinctly different from the acceleration 
region source, discussed above. 

In addition to the spectral index evolution, the fit 
results also suggest a modest decrease of the magnetic 
field and the thermal number density at the radio 
source. This behavior can be understood if we sup- 
pose, as is likely, that the flaring loop is nonuniform 
with height in such a way that higher fleld lines are 
linked with more tenuous thermal plasma (implying 
slower Coulomb loss rate) than more compact lower 
field lines. In this case the trapped fast electrons will 
live longer at the outer field lines, so the outer looptop 
regions (with lower magnetic field and smaller ther- 
mal density) will dominate the radio emission at the 
late decay phase. We note that this source remains at 
the same projected position during the entire decay 
phase as seen from the OVSA imaging. The derived 
thermal number density in the radio source is always 
below the thermal number density derived from the 
emission measure determined from the RHESSI X- 
ray spectral fit, which is consistent with significant 
spatial displacement between these two sources. 

The described fit results are highly stable vs. vari- 
ous modifications of the fitting process including vary- 
ing weights of the data points and initial values of 
the fitting parameters used as input. In contrast, 
the fit is not so stable during the second flare peak 
around 16:24 UT, which is a direct consequence of 
broader observed radio spectra, indicative, perhaps, 
of a stronger role of source inhomogeneity at this sec- 
ond peak making the spectral flt not unique. Never- 
theless, the overall range of the derived parameters is 
consistent with that for the first flare peak, although 
no direct contribution from the acceleration region is 
detectable here. 

4. DISCUSSION AND CONCLUSIONS 
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We have demonstrated that the combination of X- 
ray and microwave imaging, spectroscopy, and timing 
may allow a firm detection of the solar flare acceler- 
ation site even when the thermal SXR emission and 
microwave contribution from a magnetically trapped 
electron population are strong. For the 2002 April 11 
event discussed here, the direct radio emission from 
the acceleration region dominates the radio spectrum 
during the early, impulsive phase of the radio emis- 
sion at 2.6-4.2 GHz. 

For a better understanding of the acceleration re- 
gion and its relationship with other emission compo- 
nents it would be highly desirable to somehow con- 
strain the 3D geometry and magnetic connectivity 
within the flaring region. The current state-of-the-art 
of magnetic field reconstruction implies using a non- 
linear force-free field (NLFFF) extrapolation of vec- 
tor magnetic field measurements at the photospheric 
level to the corona. In our case, however, no vector 
magnetogram is available, but the line-of-sight data 
only from SOHO/MDI. This makes the 3D modeling 
described below non-unique, and is only offered for 
illustrative purposes. 

To build a 3D model we use our oper0 modeling 
tool, GX_Simulator, and start from the line-of-sight 
SOHO/MDI photospheric magnetic measurements at 
the region of interest. We then have the choice to ap- 
ply either potential or linear (constant-a) force-free 
field (LFFF) extrapolations to build a 3D magnetic 
data cube with which we can plot the field lines and 
form flux tube around some of them. Our tests with 
different values of force-free parameter a show that it 
is not possible to build a magnetic structure consis- 
tent with locations of the X-ray and radio sources for 
either the potential field {a = 0) or for most LFFF 
extrapolations. However, we succeeded to find the 
required connectivity for a ~ —5.5 • 10"^'^ cm~^. We 
use the corresponding magnetic data cube for our 3D 
illustrative model, summarized in Figure [TOl The left 
column of panels shows, respectively, the perspective, 
side, and top views of the magnetic structure in the 
flaring region. The locations of the acceleration re- 
gion (star symbol) and trapping source (circle) are 
shown in relation to this structure. 

It is interesting that the trapping source matches 
perfectly the top of the flaring loop as expected, while 
the acceleration region source is displaced compared 
with the loop top toward a region where many of 
the magnetic field lines are open. Moreover, inspec- 
tion of neighboring field lines implies that formation 
of a cusp nearby to the acceleration region is likely. 
Speculating further, one might conclude that the par- 
ticle acceleration site coincides with or is adjacent 
to the magnetic reconnection site. The presence of 
open magnetic field lines is also supported by lower- 
frequency radio observations including WIND detec- 
tion of type HI bursts indicating escaping energetic 
electrons. 

To complete the 3D modeling we populated the flar- 
ing loop (i.e., the region of the closed field lines shown 
by green with the central field line shown by red) with 
appropriate nonuniform spatial distributions of the 

* [www. Imsal . com/solarsof t/ssw_packages_lnf o .htmll 







Fig. 10. — 3D model of the flaring region based on LFFF 
{a Rj —5.5 ■ 10"^'' cm~^) extrapolation of the photospheric 
SOHO/MDI magnetogram. (a-c) Extrapolated 3D magnetic 
data cube visualized by two magnetic flux tubes (central fleld 
lines are red): one consists mainly of closed field lines (green) 
with a few outer open field lines (yellow), while the other one 
consists of open fleld lines only. The locations of the accel- 
eration region (star symbol) and trapping source (circle) are 
shown, (a) perspective view with the line-of-sight magne- 
togram shown at the bottom boundary; the dark red structure 
represents the model spatial distribution of the thermal plasma. 
(6) side view of the structure; step-like grey stripes represent 
the assumed hydrostatic plasma distribution in the corona out- 
side the flaring loop, (c) top view of the structure (thermal 
distribution not shown); the greyscalc background shows the 
HXR image taken at 16:20:00-16:21:00 UT. (d) top view of the 
thermal plasma distribution, (e) nonthermal electron distribu- 
tion, and (/) a combination of simulated and observed images. 
The same HXR image is used as a background in all these 
panels. The contours at the right bottom panel are for 30%, 
70%, and 90% of the corresponding peak value. Centroid of the 
acceleration region (the OVSA image at 16:20:03 UT), which 
defines the location of the star symbols in (a-c), is separately 
shown by the green triangle. 

thermal plasma (red in panels a and d) and nonther- 
mal electrons (rainbow colors in panel e). We then 
use this 3D model to compute the corresponding radio 
and X-ray images by solving the corresponding radi- 
ation transfer equations numerically. Results of this 
modeling in the form of simulated images at 22.8 keV 
and 3.5 GHz are overplotted on top of the observed 
radio and X-ray images. One can clearly see that the 
locations of the simulated images are consistent with 
each other and with observed locations of the HXR 
source and late-stage (trapped) radio source, while 
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the simulated radio source from the trapped compo- 
nent is clearly displaced compared with radio location 
of the acceleration region. 

The global parameters of the acceleration region 
are B ^ 120 G and ~ 6 • 10^'' cm~^, which arc, 
respectively, two times and ten times larger than re- 
ported in the earlier col d, tenuous (2002 July 30) flare 
(|Fleishman et al.ll2011|) . The residence time of the 
fast electrons at the acceleration region is 3 s, 
which is comparable to that in the cold, tenuous 
flare, and is much longer than the free-streaming 
time through th e acce leration region. As noted in 
[Fleishman et al.l ()201lD . this favors diffusive electron 
transport due to their scattering by turbulent waves 
and, thus, a stochastic acceleration mechanism. The 
available data is, unfortunately, insufficient to firmly 
specify the ver sion of stochast ic acceleration mech- 
anism (see e.g. lPetrosiaiill2012L for a recent review) 
operating in the event; however, it does favor those 
models predicting a roughly energy-independent dif- 
fusion time at the sour ce, like in the cold flare event 
(jPleishman et al.ll2011[ ) 

Properties of the accelerated electron components 
are somewhat different from those in the 2002 July 
30 cold flare. Pirstly, in the 2002 April 11 event 
the accelerated electron spectrum is noticeably softer 
{S « 5) than in the cold flare {S w 3.5). Secondly, in 
the cold flare the accelerated electrons arc detected 
at the energies above 6 keV, while in the April 11 
event they are only seen above ~ 20 keV; lower-energy 
X-ray emission is dominated by the thermal back- 
ground. Thirdly, the acceleration efficiency is differ- 
ent: in the cold flare almost all available electrons 
were accelerated, while in the April 11 event even the 
peak instantaneous number density of the fast elec- 
trons (rir ~ 2 • 10* cm~^) does not exceed 10% of the 
thermal electron density. Fourthly, in the April 11 
event we clearly see a spectral evolution indicative of 
the growth of a power-law tail (-Emax increases with 
time at the acceleration stage), whereas no spectral 
evolution was detected in the cold flare event, which 
implies a nearly instantaneous growth of the power- 
law tail. And finally, the released fiare energy divides 
between the thermal and nonthermal components in 
remarkably different proportions in these two events. 

Let us discuss from whence all these differences 
could originate. We have already concluded that the 
bulk acceleration mechanism is likely to be a stochas- 
tic/Permi process with a relatively long residence time 
of the electrons controlled by their spatial diffusion on 
the turbulent magnetic irregularities at the accelera- 
tion region. Por a diffusive Permi acceleration process 
the shape of the particle energy spectrum depends 
primarily on the ratio of two key parameters — the 
acceleration rate (this is the time needed to es- 
tablish the nonthermal particle spectrum, not to be 
interpreted as a duration of the acceleration process) 
and the residence/diffusion time of the electrons at 
the acceleration region in such a way that the larger 
the Ta/r^ ratio the steeper ( softer) the accelerated 
electr on spectrum, (see, e.g., [Hamilton fc PetrosianI 
|1992[) . The residence times, ~ 3 s, are comparable 
in the two events under comparison; the acceleration 
rates are, however, different. Indeed, the acceleration 



time Ta can be roughly estimated as the time needed 
for the power-law tail to grow, which is clearly shorter 
than the residence time, Tq < 3 s, in the cold flare (re- 
call, no spectral evolution was noted), while longer, 
Tq > 3 s, in the April 11 event (-Emax increases with 
time). Thus, for other conditions being equal, the 
accelerated electron energy spectral index must be 
larger in the April 11 event in agreement with obser- 
vations. 

The acceleration efficiency and energy balance in 
the flare depend, in addition to the acceleration mech- 
anism itself, on the process of electron extraction from 
the thermal pool and their injection into the main ac- 
celeration process. In the cold flare almost all avail- 
able thermal electrons were injected and accelerated, 
although their consequent energy losses were insuf- 
ficient to significantly heat the thermal plasma. In 
contrast, in the April 11 event, only a relatively mi- 
nor fraction of the thermal electrons were accelerated, 
making the collisional heating of the thermal plasma 
even less efficient than in the cold flare case (given 
that other relevant physical parameters are similar in 
these two cases). Thus, the presence of a very hot fiar- 
ing SXR plasma with T ~ 20 MK (which is present 
even before the fiare impulsive phase) requires an- 
other heating mechanism distinct from the collisional 
plasma heating by accelerated electrons. This con- 
clusion is further supported by the spatial displace- 
ment between the thermal SXR source and nonther- 
mal coronal HXR and microwave sources. 

Although the available data are insufficient to 
firmly identify the flare energization process in either 
event, or the mechanism of energy division between 
the thermal and nonthermal components, we can con- 
clude that this process does show some resemblance 
to that controlling the balance between Joule heating 
and runaway electrons in a DC electric field. Indeed, 
suppose that there is a relatively weak sub-Dreicer 
electric field directed along the flaring loop magnetic 
field. This electric field will initiate an electric cur- 
rent, which will lose its energy through Joule heat- 
ing, while the fraction of the runaway electrons avail- 
able for further stochastic acceleration will be rela- 
tively minor. In the case of a stronger electric field, 
e.g., comparable to the Dreicer field, the fraction of 
the runaway electrons becomes large, while the Joule 
heating is reduced so the plasma heating is modest. 
Even though it is a long way from these speculations 
to even a qualitative model, the analysis performed 
favors a flare picture in which electrons are flrst ex- 
tracted from the thermal pool by a DC electric field 
(of yet unspecified origin) and then stochastically ac- 
celerated to form a powcr-law-likc energy distribu- 
tion. Therefore, a stochastic acceleration mechanism 
naturally containing a DC electric field is called for. 

We have shown that observations of radio emission 
directly from the acceleration site provide important 
constraints on the acceleration mechanism in solar 
fiarcs. Despite the great differences between this fiare 
and the cold fiare, a similar acceleration mechanism, 
although operating in a somewhat different parameter 
regime, seems to be called for. Future radio spectral 
imaging observations that can better separate the ac- 
celeration site from the sites of trapping and precipi- 
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tation arc needed to investigate the flare acceleration 
niechanisni(s) in more detail. 
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Fig. 11. — Radio to HXR timing with Phocnix-2 and OVSA at different frequencies. Black thin curve, same throughout all 
the panels: HXR (20—40 keV) light curve with 2 s resolution normalized to 1; blue thick curves: resamplcd (to the RHESSI 2 s 
resolution) Phoenix-2 light curves at different frequencies indicated at the panels; green curves: spline interpolated OVSA light 
curves from original 4 s resolution to the RHESSI 2 s resolution at different frequencies indicated at the panels. The insets: 
corresponding Radio-to-HXR lag correlations in which the time delay is printed by either blue (for Phoenix-2 data) or green (for 
OVSA data). 
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Fig. 12. — Radio to HXR timing; continued. 
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Fig. 13. — Radio to HXR timing; continued. 
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Fig. 14. — Radio to HXR timing; continued. 
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Fig. 15. — Radio to HXR timing; continued. 
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Fig. 16. — Radio spectral fit of the OVSA spectra; Part II. 




Fig. 17. — Radio spectral fit of the OVSA spectra; Part III. 
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Fig. 19. — Radio spectral fit of the OVSA spectra; Part V. 



